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Abstract 

The applied method of the amplitude envelopes give us the possibility 
to describe a new class of amplitude equations governing the propagation 
of optical pulses in media with dispersion, dispersionless media and vac- 
uum. We normalized these amplitude equations and obtained five dimen- 
sionless parameters that determine different linear and nonlinear regimes 
of the optical pulses evolution. Surprisingly, in linear regime the nor- 
malized amplitude equations for media with dispersion and the amplitude 
equations in dispersionless media and vacuum are equal with precise - ma- 
terial constants. The equations were solved for pulses with low intensity 
(linear regime), using the method of Fourier transforms. One unexpected 
new result was the relative stability of light bullets and light disks and 
the significant reduction of their diffraction enlargement. It is important 
to emphasize here the case of light disks, which turn out to be practically 
diffractionless over distances of more than hundred kilometers. 

PACS 42.81.Dp;05.45.Yv;42.65.Tg 

1 Introduction 

Laser physics experimentalist have reached the point of producing without par- 
ticular difficulties few picosecond or femtosecond (fsec) optical pulses with equal 
duration in the x, y and z directions, the so called Light Bullets (LB), and fs 
optical pulses with relatively large transverse and small longitudinal size (Light 
Disks or LD). The evolution of so generated LB and LD in linear or nonlinear 
regime is quite different from the propagation of light filaments and light beams 
and has drawn the researchers' attention with their unexpected dynamics be- 
havior. For example, in pQ self-channeling of femtosecond pulses in air below 
the nonlinear collapse threshold, i.e. in linear regime, was observed. This is 
in contradiction with the well known paraxial approximation, which gives good 
agreement between theory and experiment for optical beams or long pulses with 
picosecond time duration. The main characteristic of the paraxial approxima- 
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tion is the use of a transverse Laplacian and the neglecting of the longitudinal 
second derivative of the amplitude function in space. This leads to Fresnel's 
type of diffraction for optical beams and long pulses in linear regime. Vari- 
ous different nonlinear generalizations of the paraxial approximation have been 
suggested based on the so-called spatio-temporal evolution equation [21 El E] ■ 
We suggest here another kind of nonparaxial generalization of the amplitude 
equations in dependence of the characteristics of the media (s(lo) and 712), the 
spectral region (ko and u)q), and the shape of the pulse. To provide a correct 
analysis, we need to introduce the different spatial size of the optical pulses in 
the normalized equations. The initial longitudinal spatial size is determined 
from the initial time duration by the relations zq = v g to = cto/n g , where v g is 
the group velocity and n g is the group velocity refractive index. For example, 
the typical transverse size of laser pulses is evaluated from rj_ — 1 — 5 mm to 
rj. = 100 fim 7 while for the pulses longitudinal size or the spatial period of mod- 
ulated waves, there are mainly three possibilities: a. Pulses from nanoseconds 
to 30 — 40 picoseconds. For such pulses zq >> r± and their shape is closer to 
light filaments (LF). b. Pulses from few picoseconds to 300-400 fs. In this case 
zq ~ r± and the pulses resemble LB. When the transverse and longitudinal parts 
are approximately equal (case b), it is not possible to ignore the second deriva- 
tive along the z direction in the amplitude equation, c. Pulses from 200 — 300 
fs to 10— 15 fs and less. In this case the transverse size exceeds the longitudinal 
one Zq << rj_ and the pulses initial shape can be approximated with light disks. 
The propagation of such a type of short optical pulses in nonlinear media must 
be characterized not only by a different kind of dimensions in the x, y and z 
directions, but also by a non-stationary optical linear and nonlinear response 
of the media, which lead to linear and nonlinear dispersion of different orders. 
For example, the diffraction length Zdiff = fco r i_ of LB and LD is of the order 
of the dispersion length Zdi S = ^o/fc" f° r s °lids materials and liquids in the UV 
region. This is why we should expect that the equations governing the propaga- 
tion of such types of pulses, LB and LD, would be different from the well known 
paraxial approximated equations, which describe the evolution of optical beams 
[SI El E| and nanosecond or hundred picosecond pulses • The inclusion of 
a longitudinal second derivative in space and of second-order dispersion terms, 
which are not small for LB and LD, changes the main result dramatically: even 
in the case of a linear regime of propagation, the LB generated are relatively 
stable and the LD are extremely stable in comparison with long pulses and light 
beams. 

2 From the Maxwell's equations of a source- 
free, dispersive, nonlinear Kerr type medium 
to the amplitude equation 

The propagation of ultra-short laser pulses in isotropic media, can be charac- 
terized with the following dependence of the polarization of first Pu n and third 
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P n i order on the electrical field E: 



(S{t - r) + 4ttx (1) (t - rfj E (r, r)dr = 

(< - r) E (r,x,y,z)dT, (1) 



t t— Ti t — ri— t 2 



Pi? = 3* / / / x (3) («-n,t-r 2l f-T 3 ) 



— oo — oo -co 



x (£(n,r) • B*(r 2 ,r)) E^^dr^drs, (2) 

where x^ and £ are the linear electric susceptibility and the dielectric constant, 
X*- 3 ' is the nonlinear susceptibility of third order, and we denote r = (x,y,z). 
We use this expression of the nonlinear polarization J2J), as we will investigate 
only linearly or only circularly polarized light. The Maxwell's equations in this 
case become: 



IdB 

Vx E =--_, (3) 



V • D = 0, (5) 
V • B = V • H = 0, (6) 

B = H, D = P Un +P nl , (7) 

where E and H are the electric and magnetic fields intensity , D and B are the 
electric and magnetic inductions. We should point out here that these equations 
are also valid when the time duration of the optical pulses to is greater than 
the characteristic response time of the media To (to >> To), and when the time 
duration of the pulses is of the order of time response of the media (to < To). 
Taking the curl of equation J3J and using 10} and Q), we obtain: 



^(^^)-^ = -~, (8) 
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where A = V 2 is the Laplace operator. 

Equation |JSJ is derived without using the third Maxwell's equation. Using 
equation JSJ and the expression for the linear and nonlinear polarizations Q 
and (O, we can estimate the second term in equation JSJ for arbitrary localized 
vector function of the electrical field. It is not difficult to show that for localized 
functions in nonlinear media with and without dispersion V--E = and equation 
JSJ) we can write: 

We will now express the electrical field in linear and nonlinear polarization on 
the right-hand side of with it's Fourier integral: 



+00 

E(r,t)= J E(r,u)exp(-iui)dw, (10) 
—00 

We thus obtain: 
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AE = —-^ \ J e(t — t)E (r, u) exp (—icjT)dwdT 



t t—Tt t—Tl —T-2 

' c? dt 2 



-OO —OO —OO 



X W {t - n,t - T 2 ,t - t 3 ) E(r,w) E(r,tu) 

x exp (— i (w(ri - t 2 + T 3 )))dTidT 2 dT 3 . (11) 

The causality principle imposes the following conditions on the response func- 
tions: 



e(i - r) = 0; X {3) (* - n, t - r 2 ,t - r 3 ) = 0, 

t-r<0; t-n< 0; i = 1,2,3. (12) 

This is why we can extend the upper integral boundary to infinity and use 
the standard Fourier transform jl()| : 



t +00 
J e(r — t) exp (iuJT)dT = J e(r — t) exp (iuJT)dT (13) 

— 00 — CO 

t t—Tl t—Tl —T2 

J J J X (3) (t - Tl, t - T 2 , t - T 3 ) dT X dT 2 dT 3 = 
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III X(3) ( i-Tl ' t-T2 ' < ~ T3 ) drirfT2dT3 - 



(14) 
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The spectral representation of the linear optical susceptibility £q(lu) is con- 
nected to the non-stationary optical response function by the following Fourier 
transform: 



i(u>) exp (— icot) 



:(t — t) exp (~iu!T)dT. 



(15) 



The expression for the spectral representation of the non-stationary nonlinear 
optical susceptibility is similar : 



-(-00 +00 +00 



x (3) Hex P Hwt) = - J J J x {3) (t-n,t-T 2 ,t-r 3 ) 



-00 —00 —00 



x exp (—i (w(n - r 2 + T 3 )))dTidT 2 dT 3 . 
Thus, after brief calculations, equation Ijlll) can be represented as 
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AE = 



uj 2 e{uj) 



E (r, ui) exp {—iujt)duj 



u 2 X {3) (w) 



E(r,uj) E (r,to) exp (—i (ujt))duj. 
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We now define the square of the linear k 2 and the nonlinear fc 2 ; wave vectors as 
well as the nonlinear refractive index n 2 with the expressions: 



,2 u 2 e(lu) 



(18) 



3ttw 2 x (3) H 



k 2 n 2 , 



(19) 



where 



n 2 (w) 



3ttx (3) (w) 



(20) 



In terms of these quantities, the equation Ijl7|l can be expressed by: 
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AE 



k 2 (lj)E (r, cj) exp {—iujt)duj 



k 2 (ui)n2(uj) E(r 7 uj) E (r,u) exp (—i (u>i))duj. 



(21) 



Let us introduce here the amplitude function A(r, t) for the electrical field 
E{r,t): 



E (r, t) = A (r, t) exp (i(koz — u>ot)) + c.c, 
the Fourier transform of the amplitude function A(r, lu — loq): 



(22) 



+ 00 



A (r, £) = / A (r, cj — cjq) exp (— i(u; — uJo)i)duj, 



(23) 



and the following relation between the Fourier transform of electrical field and 
the Fourier transform of the amplitude function: 



E (r, uj) exp(— iuit) = exp (— i {k^z — wot)) A (r, cj — u>q) exp — wo)i), (24) 

where wo and fco are the optical frequency and the wave number of the optical 
field, respectively. Since we investigate optical pulses, we assume that the am- 
plitude function and its Fourier expression are time - and frequency-localized. 
Substituting (|22Jl . l|23Jl and l|24Jl into equation Ij21|l we finally obtain the following 
nonlinear integro-differential amplitude equation: 



AA(r, t) + 2ik ^^ - k 2 A(r, t) = 
k 2 (ui) ( 1 + ri2(w) A (r, ui — cjq) I A (r, u — luq) exp (— i{ui — u)o)t)du (25) 



The equation H25[) was derived with only one restriction, namely, that the 
amplitude function and its Fourier expression be localized functions. This is 
why, if we know the analytical expression of k 2 (u>) and n2(w), the Fourier integral 
on the right- hand side of (|25|1 is a finite integral away from resonances. In this 
way we can also investigate optical pulses with time duration io of the order of 
the optical period To = 27r /u)q. Generally speaking, using the nonlinear integro- 
differential amplitude equation (|25(l we can also investigate wave packets with 
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time duration of the order of the optical period, as well as wave packets with a 
large number of harmonics under the pulse. The nonlinear integro-diffcrential 
amplitude equation (|25|l can be written as a nonlinear differential equation for 

the Fourier retransform of the amplitude function A, after we apply the time 
Fourier transformation 1)23)1 to the left-hand side of ()25l) : 



AA (r, uj — ujq) + 2iko 



OA (r, oj — ujq) 
dz 

+ ((l + n 2 (u) l(r,oj-u; ) k 2 (uj) - Jc%(uj )\ A(r,u- uj ) = 0. (26) 

We note here the well-known fact that the Fourier component of the ampli- 
tude function in equation (|26|l depends on the spectral difference Auj = u> — ujq, 
rather than on the frequency, as is the case for the electrical field. 



3 From amplitude equation to the slowly vary- 
ing equation of amplitudes (SVEA) 

The equation l|25|) is obtained without imposing any restrictions on the square of 
the linear k 2 (uj) and nonlinear = k 2 (uj)ri2(uj) wave vectors. To obtain SVEA, 
we will restrict our investigation to the cases when it is possible to approximate 
k 2 and k 2 ^ as a power series with respect to the frequency difference ui — luq as: 



k 2 (w) 



^ 2 £o M ,2, , , ^ (fc 2 (Up)) 
2 = k (wo) H t. (uj - UJq) 



duj 

1 d 2 (k 2 (ujq)) 

2 duj 2 



(UJ - UJof + 



(27) 



& (-) - - ^ M + "-^^ (uj - .o) + ... (28) 

To obtain SVEA in second approximation to the linear dispersion and in first 
approximation to the nonlinear dispersion, we must cut off these series to the 
second derivative term for the linear wave vector and to the first derivative 
term for the nonlinear wave vector. This is possible only if the series ()27f) and 
H28[) are strongly convergent. Than, the main value in the Fourier integrals 
in equation (|25(l yields the first and second derivative terms in 1)2 7[l. and the 
zero and first derivative terms in 1)28(1 . The first term in ()27)l cancels the last 
term on the left-hand side of the equation 1)25)1 . The convergence of the series 
1)271) and 1)281) for spectrally limited pulses propagating in the transparency UV 
and optical regions of solids materials, liquids and gases, depends mainly on 
the number of harmonics under the pulses For wave packets with more 
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than 10 harmonics under the envelope, the series (|27[) is strongly convergent, 
and the third derivative term (third order of dispersion) is smaller than the 
second derivative term (second order of dispersion) by three to four orders of 
magnitude. In this case we can cut off the series to the second derivative term 
in 1)27(1 . as the higher terms in the series contribute very little to the Fourier 
integral in equation (|25[1 . In the case of 3-7 harmonics under the pulse, the 
series l|27|) is weakly convergent and we must take into account the dispersion 
terms of higher orders of as small parameters . In the case of wave packets with 
only one or two harmonics under the pulse the series H27|l is divergent. This is 
the reason why the SVEA does not govern the dynamics of wave packets with 
time duration of the order of the optical period. As we pointed out above, this 
restriction is not valid for the amplitude equation 1251) and we can use it to 
investigate such kind of pulses. Substituting the series (|27|l and (|28|l in l|25[l 
and bearing in mind the expressions for the time-derivative of the amplitude 
function, the SVEA of second order with respect to the linear dispersion and 
first order with respect to the nonlinear dispersion is expressed in the following 
in next form: 



AA + 2ik 



dA 



0" 



2ik k 



(k k" + k r 



d 2 A 



k 2 

ft 0nZ 



-4 



A 2ikn n ik n 



Onl^Onr 



,dA 

_ 2 

A 

~dt 



A 



(29) 



We will now define other important constants connected with the wave packets 
carrier frequency: linear wave vector ko = k(uo) = ajo^/e(uio)/c linear refractive 
index ti(ujo) = y/s(u)o), nonlinear refractive index n 2 (uo) = 3irx i ul o)/ £ ( ul o), 
group velocity: 



«(«*>) = h = ] TFT". (*» 



nonlinear addition to the group velocity fc^: 



2 V e dui 

2 



2k n 2 



2 dn 2 



(31) 



and dispersion of the group velocity fc"(w ) = d 2 k/doj 2 =UJg . All these quantities 
allow a direct physical interpretation and we will, therefore, rewrite equation 
l|29|) in a form consistent with these constants: 



dA 

~at 



dA 

dz 



n 2 



kov dn 2 



A 



A 



dt 
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— AA-- fk" + —) ^A + b^l A 2 A (32) 

This equation can be considered to be SVEA of second approximation with 
respect to the linear dispersion and of first approximation, to the nonlinear 
dispersion (nonlinear addition to the group velocity). It includes the effects 
of translation in z direction with group velocity v, self-steepening, diffraction, 
dispersion of second order and self-action terms. 

4 Propagation of ultrashort optical pulses in vac- 
uum and dispersionless media 

The theory of light envelopes is not restricted only to the cases of non-stationary 
optical (and magnetic) response. Even in vacuum, where e = 1 and P n i — 
we can write an amplitude equation by applying solutions of the kind (|22|l to 
the wave equation @. We thus obtain the following linear equation for the 
amplitude envelope of the electrical field: 

/ dl dl\ c - 1 d 2 A 
~\^ +C ^)-2Tv~2k^ C ^- 

The Vacuum Linear Amplitude Equation (VLAE) l|33l) is obtained directly from 
the wave equation without any restrictions by using the series of the square of 
the wave vector, as is the case of a media with dispersion, where it is required 
that the series Q27JI be strongly convergent. This is why equation (|33|l describes 
both amplitudes with many harmonics under the pulse, and amplitudes with 
only one or a few harmonics under the envelope. It is obvious that the envelope 
A in equation l|33|) will propagate with the speed of light in vacuum c. The 
equation (|33|) is valid also for transparent media with stationary optical response 
e = const. In this case, the propagating constant will be v — c/^/eJl. 

5 SVEA and VLAE in a normalized form 

Starting from Maxwell's equations for media with non-stationary linear and 
nonlinear response, we obtained an amplitude equation and a SVEA using only 
two restrictions, which are physically acceptable for ultra-short pulses. Having 
adopted the first restriction, namely, investigation of localized in time and space 
amplitude functions only, we introduced the amplitude equation l)25|) . Following 
the second restriction, i.e., limiting ourselves with the case of a large number of 
harmonics under the localized envelopes, we obtained the SVEA (|32"|) . As it was 
pointed out in the previous section, these restrictions do not affect the VLAE 
The next step is writing SVEA and VLAE in dimensionless 
variables and estimating the influence of the different differential terms. In this 
case, the coefficients in front of the differential operators in (|32|l and H33JI will 
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be numbers of different orders, depending on the media n and n 2 , the spectral 
region of propagation ko and ujq, the field intensity |A)| 2 , and the initial shape of 
the pulses, namely, optical filament rj_ << zq, optical bullets r± ss zq or optical 
disks rj_ >> zq. With r± we denote here the initial transverse dimension, "the 
spot" of the pulse, and with zq we denote the initial longitudinal dimension, 
which is simply the spatial analog of the initial time duration to, determined by 
the relation zq = vto or zq = cio in the vacuum case. The SVEA (US2[| and VLAE 
l|33l) are written in a Cartesian laboratory coordinate system. To investigate the 
dynamics of optical pulses on long distances, it is convenient to rewrite these 
equations in a Galilean coordinate system, where the new reference frame moves 
with the group velocity, t' — t; z' = z — vt for equation l|32l) : 



dA 



n 2 



kov dn 2 



\ 



A 



A I d 



A 



A 



dv 



dz' 



v 3 kl d 2 A v 



2 dz r 



k" 



d 2 A d 2 A 
~W 2 ~ 2v dt'dz' 



2k 
n 2 k v 



A±A - 



A 



2 A, (34) 



and with velocity of light for equation (|33|l . t' = t: z' = z — ct: 



~ l dF~ 2k~o ~~ 2k c dt' 2 + k Q dt'dz" ( ' 

and with Aj_ = + -§jj2 we denote the transverse Laplacian. We define the 
following dimensionless variables, connected with the initial amplitude and with 
the spatial and temporal dimensions of the pulses with the relations: 



A = A A" ; x = r^x"- y = r ± y"; z' = z z"; t' = t t"; z = z z v ; t = tot". (36) 

We substitute these variables in 13211 . (J33J, l|33l) and l|35l) and, making use of the 
expressions for the diffraction Zdif = k$r\ and dispersion Zdisp = to/k" lengths, 
we obtain the next five dimensionless parameters in front of the differential 
terms in the equations i|3"H J33J) and (|33j) : 



a = k a z Q ; S 2 = /? = 7 = k 2 rln 2 \A \ 2 ; 

**0 %disp 

7l = 2fc ron 2 |A,| 2 (n 2 + ^|^). (37) 

In the new dimensionless variables and constants, the equations l|32l) . I|34l) . 13311 
and (|3"5|) can be represented as (the seconds of dimensionless variables are omit- 
ted for simplicity): 
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Case a. S VEA (KS2I) in a laboratory frame (" Laboratory" ) 



2iaS 2 



dA dA 
dt dz 



•7i- 



.4 



.4 



\ 



dt 



d 2 A d 2 A 

W 2 ~ ~w 



(38) 







d 2 A 
~b~t 2 



A 



2 _ 

A. 



Case b. SVEA H34() in a frame moving with the group velocity: 



2iaS 2 



dA 



A 



A 



dt' 



A 



A 



dz 1 



(13 + S 2 ) 



d 2 A 
W 2 



A ± A-(3 



, d 2 A 
'dt'dz' 



d 2 A 
dz 12 



(39) 



.4 



A, 



Case c. VLAE (13311 in a laboratory frame: 
Case d. VLAE l|35|) in a Galilean frame: 



d 2 A d 2 A 



dz 2 



dt 2 



(40) 



2dA_ A x s Jd 2 A 



2ia5 z — = Aj_A 
dt' 



d 2 A 



dt' 2 dt'dz 



(41) 



We obtain equal dimensionless constants in front of the differential terms in 
both the "Laboratory" and "Galilean" frames. This gives us the possibility to 
investigate and estimate simultaneously the different terms in the normalized 
equations l|3*%)l. fflfy . and (f!T|l . We will now discuss these constants in 

more detail, as they play a significant role in determining the different pulse 
propagation regimes. 

- The first constant a = kgzo — 2nzo/Xo determines with precision 2ir the 
"number of harmonics" on a FWHM level of the pulses. Since we use the slowly 
varying amplitude approximation, a is always a large number (a >> 1). 

- The second constant S 2 = r 2 ^ / z 2 determines the relation between the initial 
transverse and longitudinal size of the optical pulses. This parameter distin- 
guishes the case of light filaments (LF) 8 2 = r 2 ^/ z 2 « 1 from the case of LB 
5 2 = t 2 ± /zq = 1 and the case of light disks LD S 2 = r\j ' z 2 >> 1. For light 
filaments 5 2 << 1 and we can neglect the differential terms with coefficient 
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S 2 . It is not difficult to see that in this case the SVEA (gHJ, ® and VLAE 
(|39() . (|40|) can be transformed to the standard paraxial approximation of the 
linear and nonlinear optics. If we set the possible values of the optical pulses 
transverse dimensions at 3 — 4mm > r± > 100/im, we can directly obtain above 
distinction in dependence of time duration of the pulses. For light pulses with 
time duration ns > t > 40 — 50ps we obtain S 2 — r\j ' z 2 « 1 and we are in 
regime of LF and paraxial approximation. In case of light pulses with duration 
to w 3 — 4ps up to 500 — 600/s it is possible to reach S 2 = t\/Zq = 1 and we 
are in regime of LB. For pulses in the time range 300/s — 30/s we can prepare 
the initial shape of the pulses to satisfy the relation S 2 = r\j ' z 2 > > 1 and thus 
reach the LD regime. It is important to note here that the wave packets in visi- 
ble and UV ranges with time duration to > 30/s admit more than 10-15 optical 
harmonics under the pulse, so that we are still in SVEA approximation. In the 
last two cases (LB and LD) the differential terms with 5 2 can't be ignored and 
the equations (|3*S)l . , jSHJl and l|4T)| governed the propagation of pulses with 
initial form of LB and LD are quite different than paraxial approximation. 

-The third parameter is /3 = fco r i/ Zdis-, where z<u s = to/k", determines the 
relation between diffraction and dispersion length. The dispersion parameter 
fc" in the visible and UV transparency region of dielectrics has values from 
k" ~ 10 -31 s 2 1 cm for gases and metal vapors up to k" ~ 10~ 26 s 2 /cm for solid 
materials. It is convenient to express this parameter using the product of the 
second constant S 2 and the parameter (3\ — kov 2 k" by the relation (3 = /3iS 2 . 
For the typical values of the dispersion fc" in the visible and UV region listed 
above, the dimensionlcss parameter /3i is very small (f3i << 1) in the visible 
part of the spectrum, while for optical pulses propagating in the UV region in 
solids and liquids it may reach (3\ oc 1. The parameter (3\ can be also negative 
and may reach (3\ ~ — 1 near electronic resonances and also near the Langmuir 
frequency in electronic plasmas ■ It is no difficult to see that only in the case 
/?i ~ — 1 and 5 2 oc 1 can we obtain the 3D+1 nonlinear Schrdinger equation 
from SVEA Qg) and © 

- The fourth and fifth constants 7 = k^r^ni \ Ao\ 2 /2 and a5 2 ^i — a6 2 n 2 \Ao\ /2 
are correspondingly the nonlinear coefficient and the coefficient of nonlinear ad- 
dition to the group velocity (coefficient in the front of the first order nonlinear 
dispersion). It is easy to estimate that for a >> 1 and S 2 > 1, we always have 
7 >> oc8 71. In this paper we investigate optical pulses with power near the crit- 
ical threshold for self- focusing 7 = 1 and less (linear regime) 7 << 1. For these 
cases the nonlinear addition to the group velocity is very small (a<5 2 7i << 1) 
and from here to the end of this paper we will neglect the terms with the first ad- 
dition to the nonlinear dispersion. The analysis of the dimensionless constants 
performed above leads us to the following conclusion: The investigation of wave 
packets with power near that for self-focusing 7 cx 1 in the visible and UV region 
in a media with dispersion are governed by the following SVEA equations: 

Case a. Laboratory frame (" Laboratory" ) 




) 



= A ± A + & 



. 2 d 2 A 



dz 2 



S 2 (01 + 1) 



d 2 A 



dt 2 



+ 7 A A. (42) 
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Case b. Frame moving with group velocity ("Galilean"): 



2ia8^ = A ± A-(3 1 S^-5*([3 1 + l) 



d 2 A d 2 A 

W 2 " 1 



dt'dz 



7 



.4 



A. (43) 



Equations 1(42(1 and 143(1 differ from the paraxial spatio-temporal evolution equa- 
tions investigated in [21 El EI I13| by the inclusion of a second derivative along 
the z direction, mixed term and additional second derivative in time term. This 
leads to a quite different dynamics of the ultrashort pulses; we will first inves- 
tigate their propagation in linear regime, when 7 << 1. 



6 Analytical solutions of the linear SVEA and 
VLAE 

The behavior of long pulses is similar to that of optical beams, since their 
propagation is governed by a similar kind of paraxial equations, as pointed 
above. That is why we cannot expect the diffraction enlargement of long pulses 
to be different from that of optical beams. The situation with the LB and LD is 
different. Their propagation is governed by the new spatio-temporal evolution 
equations in media with non-stationary optical response - SVEA 1(42(1 and ((43(1 , 
and by VLAE 1(33(1 and ((351 in media with linear stationary optical response (or 
vacuum). In this section we will attempt to solve the equations 1(42(1 . 1(43(1 in 
linear regime (7 << 1) and will compare the solutions with the solutions of the 
linear VLAE l(33|) and l(35|l . Neglecting the small nonlinear terms in l(42|) . 1(43(1 
we arrive at: 

a. Linear SVEA in a laboratory coordinate frame: 



„. c2 OA dA\ s2 d 2 A r2/a ^d 2 A 

- 2iaS \JH + ~~dz J = x + ~dz 2 ~ ^ + 1} W (44) 

b. Linear SVEA in a Galilean coordinate frame: 



- arfg - A±A - S 2 (ft + I) (p - 2^) - ^ A 0. (45) 

For comparison we will rewrite here the corresponding linear VLAE: 
c. Linear VLAE in a laboratory frame: 



n . r2 dA dA\ - , 2 d 2 A d 2 A\ 



d. Linear VLAE in a Galilean frame: 
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The important result here is that in linear regime the normalized equations for 
media with dispersion l|44l) , i|45|) identical equal with accuracy of one dimension- 
less dispersion parameter j3\ with the equations for media without dispersion 
and vacuum I^J), (gZJ. In the equations (H), G3> and igBJl, lO there are 
only three dimcnsionless parameters, namely, a = fco?"o » L which determines 
the "number of harmonics" on FWHM level, S 2 — r^/z 2 , which determines the 
relation between the longitudinal zq and the transverse r± size of the optical 
pulse, and the dispersion parameter (3\ — kov 2 k" . As it can be expected, the 
equations for ultra-short optical pulses in vacuum and dispersionless media (1461) , 
<|47|) become identical with the equations with dispersion (|44H . (|45|l when the 
dispersion parameter (3\ « 1 and the product of 8 2 (3\ « 1. 



6.1 Linear regime of long optical pulses (Paraxial approx- 
imation) 

It is easy to see that in the case of long pulses (from ns to 50 — 60ps) S 2 « 1. 
The values of the dimensionless dispersion parameter (3\ range from a small 
number in the visible region up to (3\ oc 1 — 2 in UV transparency region of 
solids, as it was pointed above. This is why we can neglect the terms multiplied 
by the small parameter 5 2 « 1 in l|44|) . I|45|l . (|46|l and (|47|l . and also the terms 
containing the product of a small parameter with a parameter of the order of 
unity <5 2 (1 + (3\) « 1 and 5 2 (3i in ifUjl . (|4"5|) . For wave packets with a high 
number of harmonics under the pulse, the condition aS 2 w 1 is usually fulfilled. 
Thus, for long pulses or LF, we reduce these equations to the well-known scalar 
paraxial approximations: 

a. SVEA and VLAE in a laboratory frame: 

+ ^)+^ = »- <«» 

b. SVEA and VLAE in a Galilean frame: 



-2iaS 2 — +A ± A = 0. (49) 

It is clearly seen that the dynamics of long pulses propagating in the trans- 
parency region of gases and solids is governed by the same equations that 
describe the evolution of long pulses in vacuum and are similar to the opti- 
cal beam equation. If we turn our attention to the dimension variables, we 
will see that there is only one main difference between propagation of LF in 
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a media with non-stationary optical response and propagation of LF in vac- 
uum. In paraxial SVEA, the long pulses propagate with group velocity v, while 
in paraxial VLAE, the long pulses propagate with velocity of light c. Let us 
present the initial condition of the paraxial equations l|48ll , (|49fl correspondingly 
A L (x,y,z + t = 0) = A%(x,y), A G {x,y,t = 0) = A$(x,y), where A$, A§ and 
Aq are the initial amplitude functions in "Laboratory" and "Galilean" coor- 
dinates. Then the solutions of equations JISJ, 1(49(1 are well known and they 
are obtained as a convolution of two inverse Fourier transforms, the Fourier 
transform of the initial pulse (in our case Gaussian pulse) and the spectral 
presentations of the Fresnel's kernel in "Laboratory": 



A L (x,y,z + t) = F- 1 A%(k Xl k y ) 



I iaS 2 

V(Z + t) J J R 2 



exp 



2(z + t) 



exp 



-i(z + t) 



2ia- 



(50) 



[{x-x'Y + {y~y'?])y-A L {x> iy ') 



dx' dy' , 



'Galilean": 



A L (x,y,t') = F~ 1 A^(k x ,k y ) 



IF" 1 



exp 



(51) 



ia5 2 
nt' 



exp 



iaS 
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r [(x~x') 2 + (y-y') 2 } HVi) 



dx' dy'. 



As it can be seen from the solutions in the different coordinate frames, we have 
exact Fresnel's integrals and phase modulation not only in time ("Galilean") 
H51|) but also in space-time ("Laboratory" )J50J|. These considerations provide 
an insight into the real phase modulation of long optical pulses in paraxial 
approximation. The pulses obtain equal amounts of space and time phase mod- 
ulation under propagation in the media and in vacuum. This is an important 
difference from beam propagation, where there is no group velocity or velocity 
of light in vacuum in the laboratory frame and the beam undergoes only spatial 
modulation. We can therefore make the following conclusion: The propaga- 
tion of long pulses in the region of transparency of a media and in vacuum are 
equal with accuracy of the type of velocity of propagation and undergo Fresnel's 
diffraction. 



6.2 LB (from few ps to 300-400 fs) and LD (from 200 fs 
to 30-40 fs) in linear regime 

We will apply the Fourier method to solve the linear SVEA l|44|l , 145(1 which 
describe the propagation of ultrashort optical pulses in media with dispersion 
and linear VLAE (gHJ), which govern the propagation of light pulses in 
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vacuum and dispersionless media. To distinguish between the solutions of SVEA 
and VLAE we will mark the Fourier transform of the amplitude functions of 
SVEA 1331), gSJl in a Galilean frame with A G (k x ,k y ,k z ,t) = F(A(x,y,z,t)), 
and in a Laboratory frame with A^{k x , k y , k z ,t) = F(A(x,y, z,t)), and the 
Fourier transform of the amplitude functions of VLAE Ij46(l , H47(l in a Galilean 
frame with Bc{k x , k y , k z ,t) = F(A(x,y, z,t)) and in a Laboratory frame with 
BL(k x , k y , k z ,t) — F(A(x 1 y,z,t)). Applying spatial Fourier transformation of 
the components of the amplitude vector functions A, we obtain the following 
ordinary linear differential equations in k x ,k y , k z space for SVEA: 
a. "Laboratory": 



- ^S 2 ^ = - (k x 2 + k y 2 + 5 2 {k 2 - 2ak z )) A L - 6 2 (A + 1) (52) 
b. "Galilean": 



- 2iS 2 (a - + 1) k z ) = - {k x 2 + k v 2 - 6 2 ^kl) A G - S 2 (ft + 1) (53) 

and the following equations for VLAE: 
a. "Laboratory": 

- 2*a6 2 ^ = {k 2 + ky 2 + 5 2 {k 2 - 2ak z )) B L S 2 ^. (54) 
b. "Galilean": 



■2^ 2 { n. ,,,) ^ (k x 2 + k y 2 ) B G S 2 ^. m 



We will now denote the square of the sum of the wave vectors with k 2 — k x 2 + 
k y 2 +5 2 (k z 2 — 2ak z ). We look for solutions of the kind Al — AL(k x ,k y ,k z ) exp^r^i) 
and Aq = Aa(k x , k y , k z ) exp(i£lct) for the equations l|52|l . l|53|l . and for solu- 
tions of the kind Bl = B^{k Xl k yi k z ) exp(i$£t) and Bq = Bc{k x , k y , k z ) exp(i^ct) 
for the equations l|54|) and l|55|) correspondingly. The solutions exist when 
D,g, $l and <&g satisfy the following quadratic equations: 



o2 , (a-(^i + l)M n k x 2 + k y 2 ~5 2 p 1 k 2 z 

" G + 2 ]h + i WWTV) — = - (57) 



1G 



$1 + 2a<S> L -^2=0, (58) 



$a, + 2(a-fc,)$ G - *^±*g! =0. (59) 



The solutions of , JSTJ) , (JSHJl, and (J55J) are similar: 



i o a , I a 2 k 2 



nL ' = -—i ± V wn? + mTT)' (60) 



12 g-^ + l)^ / (a-Qgi + l)*;,) 2 fc 2 + fc 2 - £ 2 /?ifc 2 

= aTi V (/3 1 + D 2 + WTI) ' (61) 



$ L M = -aiVa 2 + fc 2 /« 2 , (62) 



$ G ^ = -(a-jy iyV+P/tf 2 . (63) 

Now the necessity of a parallel investigation of the propagation of optical pulses 
in media with dispersion, in dispersionless media and in vacuum becomes obvi- 
ous. Further, we will introduce here the concept of week and strong dispersion 
media depending on the value of dimcnsionlcss dispersion parameter pi . When 
Pi « 1 we have media with week dispersion. It is no difficult to calculate that 
in optical transparency region of gases, liquids and solids materials, this param- 
eter is usually is very small ( pi << 1). For media with week dispersion the 
solutions of the characteristics with dispersion i|tj(Jfl , (|61f) are identical with the 
solutions for media without dispersion and vacuum i|62|) . I|63[) correspondingly. 
In media with strong dispersion Pi can reach the values Pi ~ 1 — 3 in the UV 
transparency region of solids and liquids. In this case, the solutions for media 
with dispersion will be slightly modified with the factor (3% + 1 with respect to 
the solutions without dispersion. We consider here the regime of propagation 
far away from electronic resonances and the Langmuir frequency in electronic 
plasmas, where it is possible to obtain a strongly negative dispersion parameter 
Pi oc —1. We point here again, that in the case of LB, when f3% = —1, the am- 
plitude equations l|32|) can be transformed into the 3D+1 linear and nonlinear 
vector Schrdinger equations |14j . Generally speaking, the dispersion parameter 
Pi varies slowly from the visible to the UV transparency region of the materials 
from very small values up to Pi ~ 1 — 3; this is why it does not influence radi- 
cally the solutions and the propagation of optical pulses in linear regime. The 
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other parameters a and S 2 change significantly. For example a varies from 10 1 
to 10 3 while S 2 varies from 1(T 2 - 1(T 4 for LF to 10° for LB and 10 2 - 10 4 for 
LD. This is why in the next paragraph we will investigate more precisely the 
solutions of the equations for media with week dispersion, dispersionless media 
and vacuum 1)62)) . 1)63)) . as we expect that the solutions for media with strong 
dispersion (UV transparency region of solids and liquids) i)6(J)l . i)61|) will be only 
slightly modified by the factor fi\ + 1. As we obtained the solutions for the 
characteristics (|60|) . 1)6 ljl. 1)62)1 and 1631) . the solutions of the corresponding lin- 
ear differential equations SVEA lt55|l and VLAE J53J) in the fc-space 
become: 

a. Solution of SVEA in the k-space and a laboratory coordinate frame: 



A L = A L (k x ,k y ,k z ,t = 0) x 







' a 2 


k 2 



exp - «4t ± W -t^tt + " | ^ 1*1- (64) 



b. Solution of SVEA in the k-space and a Galilean coordinate frame: 

A G = A G (k x ,k y ,k z ,t = 0) x 



, , a-(pi + l)k z / (a - [fa + l)k z ) 2 k 2 x + k 2 ~ 5 2 ^k 2 z 

exp| - z| — pm V (^+d 2 m+D 1 '- ( } 

c. Solution of VLAE in k-space and laboratory coordinate frame: 



B L =B L (k x ,k y ,k z ,t = 0) exp I -a ± V" 2 + ^V^J * 1 • (66) 
d. Solution of VLAE in the k-space and a Galilean coordinate frame: 



B G = B G {k x ,k y ,k z ,t = 0) exp (-i ((a - k z ) ± ^a 2 + k 2 /5 2 ^j tj . (67) 

It is natural that the solutions l|66|l and (|67|l of equations (|54|) and (|55|l should 
identical be equal with accuracy to a wave number in the z direction. This 
follows from the Fourier transform of such evolution equations and leads to only 
one difference between the solutions in the real space - the motion of the pulse 
in the z direction in a "Laboratory" frame and its stationarity in a "Galilean" 
frame. We deal here only with initial functions for the amplitude envelopes 
of the electric field that are localized in space and time. Thus, the images of 
these functions after Fourier transform in the k x ,k y , k z space also are localized 
functions. The solutions of our amplitude equations in k space 1)64)1 . 1)65)1 . 1)66)1 
and 1)67)) are the product of the initial localized in k x ,k y , A: z -space functions and 
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the new spectral kernels, which are periodic (different from Fresnel's one). The 
product of a localized function and a periodic function is also a function localized 
in the kx •> ky ■ k z space. Therefore, the solutions of our amplitude equations in k 
space (IBT)|) . <|fi7|) . iKH|l and also are localized functions in this space and we 
can apply the inverse Fourier transform to obtain again localized solutions in 
the x,y,z,t space. More precisely, we use the convolution theorem to present 
our solution in the real space as a convolution of inverse Fourier transform of 
the initial pulse with the inverse Fourier transforms of the new spectral kernels: 
a. Solution of SVEA (|44H in a laboratory coordinate frame: 



A(x, y,z,t)=F 1 (A L (k x ,k v ,k z ,t = 0)) 



b. Solution of SVEA (|45|l in a Galilean coordinate frame: 

A(x,y,z',t') ^F- 1 (A G (k x ,k y ,k z ,t = 0)^) 
( . ( a-Q8i + l)fc« / (a - (ft + l)k z ) 2 ' kl + k 2 - 5^kj 

a. Solution of VLAE l|46|l in the laboratory coordinate frame: 

A(x,y,z,t) = F- 1 (l} L {k x ,k y ,k z ,t = 0)) ® 

F- 1 ^cxp (i (-a ± ^a 2 + k 2 /5 2 ^J tj ^ . (70) 

b. Solution of VLAE l|47|l in a Galilean coordinate frame: 



k 2 




A(x,y,z\t') = F- 1 (s G (fc x ,fc y ,fc z ,t = 0)) 



F (^exp I (a - k z )± \j a 2 + k 2 /S 2 j tjj , (71) 

where with F^ 1 we denote the spatial three-dimensional inverse Fourier trans- 
form and with (g> we denote the convolution symbol. The difference between 
the Fresnel's integrals (|50|l and l|51l) . describing propagation of optical beams 
in linear regime of propagation, and the new integrals (|68[) . <|t)9|) . (|7U|I and (|71|) . 
solutions of the linear evolution equations (|44l) . I|45l) . (|46ll and l|47l) is quite 
obvious. In addition in the new spectral kernels there are three dimensionless 
parameters, a, which determines the number of harmonics under the pulse with 
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accuracy 27r, S 2 , which gives the relation between the initial longitudinal and 
transverse dimensions of the optical pulse, and /3i, which accounts for the influ- 
ence of the second-order dispersion. We will focus mainly on the SVEA, and will 
fixed a to be always a large number a » 1. As it was pointed out above, the 
condition a >> 1 is not necessary for VLAE so that we can investigate pulses 
with longitudinal duration of order of the carrier wavelength in vacuum and in 
dispersionless media. To analyze more precisely the influence of the other two 
parameters, S 2 and /3i on the evolution of the initial pulse we will rewrite the 
expression for the spectral kernel l(fiUI) of the solutions (|fJ8|) of equation in 
the following form: 



tlr 1 ' 2 



a 2 I 1 2a 

± V Wi + 1) 2 + m + + k2y) + T^ kl " K+i kz ' (72) 

As a is a constant, the diffraction widening will be determined by the second 
and third term under the square root in 172(1 : 

which determines the transverse and longitudinal diffraction and dispersion 
widening of the pulses. We pointed out above that the dispersion parameter 
varies very slowly within the limits < (3\ < 10 1 , while the relations between 
the transverse and longitudinal part varies significantly I0 -4 < S 2 < — fO 4 . This 
is why we will estimate mainly the influence of the different values of S 2 and 
a 2 on the diffraction widening. We will consider investigate the following basic 
cases: 

a/Long pulses, when S 2 << 1. ft is easy to estimate from {73J, that the 
transverse enlargement k x 2 +k y 2 will dominate significantly over the longitudinal 
one k 2 as: 

sWTT) >;> 1; ATI ~ 1/2 - (74) 

A more precise mathematical analysis will give us the standard Fresnel's spectral 
kernel and paraxial approximation when (3% « 1. In the case of long pulses we 
have also aS 2 ~ I and it is seen from equations 1(52(1 and 1(53(1 that we can obtain 
a similar diffraction length {zdiff — fco^l), as in the case of an optical beam. 
The difference is only with the factor j3\ + 1 , while for the cases of media with 
week dispersion, dispersionless media and vacuum, there is exact coincidence 
between the diffraction of long pulses and diffraction of optical beams. 

b/LB: 5 2 ~ 1 . In the case of optical pulses with relatively equal transverse 
and longitudinal size, we obtain equal coefficients in front of the transverse 
k x 2 + k y 2 and longitudinal k 2 diffraction terms: 

nsT+I) " 1/2; aTT** ~ 1/2 ' <75) 
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Hence, the diffraction and dispersion transverse enlargement will be reduced 
with the factor 8 2 ((3i + 1) with respect to the diffraction of long pulses and 
Fresnel's diffraction. The LB widening will be further reduced if the square of 
the number of harmonics multiplied by Air 2 - a 2 exceeds the square of the sum 
of the normalized wave vectors (l/2)(/c 2 + k 2 + k 2 ). For example, the main part 
in the Fourier integral of one normalized localized Gaussian pulse gives values 
of the wave vectors in the range —4 < k Xy y iZ < 4 and the sum is of the order 
of (1/2) (k 2 + k 2 + k 2 z ) ~ 20 — 30. For a pulse with many harmonics under the 
envelope, a 2 can reach a 2 ~ 10 2 " 4 . Thus, the square of the number of harmonics 
a 2 under the pulse starts to prevail over the diffractive terms k in the Fourier 
integrals and Therefore, the two conditions, 8 2 ~ 1 and a 2 » 1, 

open up a unique possibility for a considerable reduction of the diffractive effects 
in the dynamics of LB and for a relatively stable linear propagation of LB as 
compared with long pulses and light beams. 

In addition, we will point here to an important asymptotical behavior of 
LB: When a 2 is small (pulses with only one two harmonics under the envelope) 
and 0i « 1 (media with week dispersion, dispersionless media and vacuum), 
the spectral kernels of the new equations tend to the asymptotical value ~ 
cxp [i{^jk 2 + k 2 + (k z + a) 2 )tj = cxp (i(\k\t)), which is actually the spectral 
kernel of the 3D wave equation. This is why for optical pulses with only one 
or two harmonics under the envelope we can expect diffraction similar to the 
typical diffraction of the 3D wave equation, with the appearance of observing 
of internal and external fronts and a significant widening of the pulse. 

c/ Light disks: This is the case when the longitudinal size Zo is mush shorter 
than the transverse size r± and S 2 >> 1. As it was indicated above, the typical 
time region for such pulses is 30 — 40/s < to < 200 — 300/s. We determine the 
lower limits of this relation from the conditions a 2 >> 1, i.e a large number 
of harmonics under the envelope. This condition still holds true for pulses in 
the visible and UV regions with time duration 30 — 40/s. The dimensionless 
parameter in front of the transverse diffraction and dispersion [k x 2 + k y 2 ^j will 
be of the order of: 

' « 1, (76) 



S 2 (Pi + 1) 

while for the coefficient in front of the longitudinal diffraction k 2 we have 



(ATT) al/2 - (77) 

We thus see that the transverse enlargement is of the order of S 2 (f3i + 1), or 
negligible as compared with LB, and smaller by a factor of about (<5 2 (/3i + l)) 2 
than in the cases of long pulses and Fresnel's diffraction. The longitudinal 
enlargement k\ is of the order of that for LB. To summarize the results of 
this section, we can we expect that the transverse diffraction and dispersion 
enlargement of LB should be smaller by a factor <5 2 (/3i + 1) than those of LF, 
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while the transverse diffraction and dispersion enlargement of LD should be 
smaller by a factor of {S 2 (f3i + l)) 2 than those of long pulses and paraxial 
approximation. Practically no transverse enlargement of LD would be observed 
over long distances, namely more than tens or hundreds of kilometers. 



7 Dynamics of LF, LB and LD in media with 
weak dispersion, dispersionless media and vac- 



In the previous section we first found the solutions of SVEA (|46|l and VLAE 
<|47l) in the k space (|64|I . (|65|) . H66|I and (|67|l . The dynamics of the solutions under 
different initial conditions in the (x,y,z,t), space can be found by applying an 
inverse Fourier transform to the solutions in the k space <|ol3|) . J7UJ), <|fi§|l and 
J7TJ. As it was pointed out before, the dispersion parameter varies very slowly 
in the limits < (3\ < 10 1 , while the relations S 2 and a vary significantly. We 
will, therefore, study the VLAE, as the solutions of SVEA for media with strong 
dispersionJHHJ and ljTU|) differ from the solutions of VLAE (|rT§|) and ifTTJl only by 
the factor 0± + 1. We should remind here that we investigate in this paper only 
the case of the region of transparency of the materials with positive dispersion. 
In the beginning we will fix the dimensionless parameter a (number of harmonics 
multiplied by 2ir) to be a large number number a = 40 and consider the following 
basic cases: dynamics of LF (S 2 << 1), light bullets (S 2 — 1 and light disks 
S 2 >> 1. We compare the evolutions of LF, LB and LD with the solution of 
a Gaussian beam in the framework of the standard paraxial equation of linear 
optics, written in a laboratory coordinate frame (|48|1 with initial condition of 
the type: 



The evolutions of the initial condition (|78|) governed by the paraxial equation 
H48|) is described by the Fresnel's integral or can be found by numerical cal- 
culation of the inverse Fourier transform of the solution in the (k x , A; y )-space 
<|5U|) . The solutions of the paraxial equation (|48|l with initial condition l|78ll on 
a normalized distance 47r/15 is illustrated in Fig.l. 

7.1 Evolution of long optical pulses (light filaments) 

The dynamics of optical pulses in media with week dispersion, dispersionless 
media and vacuum is governed by the VLAE (|46D . In the case of a long pulse 
we can express the initial conditions as follows: 



uum 




(78) 
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1=0 



1=20 %n$Q 



1=40 at 150 



Figure 1: Evolution of a Gaussian beam with initial condition A x (x, y, z = t = 

0) = exp ^— x \ governed by the 2D paraxial equation Ij48(l . The transverse 

size (the spot) grows approximately twice over the normalized distance z = 
4vr/15. 



A = A X X] a = 40; S 2 = ^ = — , 

A x (x, y,z,t = 0)= exp + ^ + . (79) 

The solution of VLAE in a laboratory coordinate frame (|70|l for an optical 
pulse whose longitudinal size is nine times as large as the transverse one l|79|l 
is calculated by using the numerical FFT technique. The result is presented 
in Fig. 2. The propagation distance is the same as in the previous case of a 
light beam 47r/15. Figure 3. illustrates the spot size (x,y plane) of the pulse. 
The analysis of the solution reveals the same evolution dynamics as in the case 
of an optical beam, with the diffraction length being approximately equal to 
the diffraction length of an optical beam Att/15. A more detailed investigation 
discloses the following natural dependence on the dimensionless parameters a 
and S 2 in the case of a long pulse: as the number of harmonics under the pulse 
a increases, together with the longitudinal size with respect to the transverse 
one, the VLAE solutions (|70() approach the paraxial equation solutions l|48() . 

7.2 Propagation of light bullets in linear regime 

The evolution of LB in media with week dispersion, dispersionless media and 
vacuum is governed by the same VLAE (1461) . The shape of the LB is symmetric 
in the x, y and z plane, so that the initial Gaussian profile can be written as: 
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1=0 



1=20 %! 150 



Figure 2: Evolution of a long Gaussian pulse with initial condition A x (x,y,z,t = 
0) = exp ^_ *W+a 2 * 2 ^ ) a = 40 . s = yg govern i ng by the VLAE ® . The xz 
plane of the pulse is presented. 



A = A x x; a = 40; 5 2 = = 1, 

/ + iv 2 + z 2 \ 
A x (x,y ) z ) t = 0)=expl 1 . (80) 

From the qualitative analysis presented in the previous section, we can expect 
a significant reduction of the widening of LB (of the order of l/<5 2 ) with respect 
to long pulses or optical beams. The solution of VLAE (|70|l for LB on a distance 
57T, calculated again by means of the FFT technique, is illustrated in Fig. 4. in 
the x, z plane. One can only see a negligible enlargement of the pulse along a 
considerable distance. Figure 4. presents the solution H70|l of the VLAE (|46|) . 
with initial conditions of the type (|80|) . If we consider a large distance and 
times z = t > 5ir, the LB will go out the grid. This is the reason why the 
next calculations of LB and LD evolution over long distances are performed 
in framework of the solutions l|71|l of the VLAE l|47|l written in a Galilean 
coordinate system under similar initial conditions: 



A = A x x: a = 40: 5 2 = = 1, 

A x (x, y, z\ t' = Q) = exp (- ^ + £ + - ) . (81) 

The LB evolution in a Galilean frame over a long normalized distance t = 15w, 
where the transverse width of the pulse is approximately doubled, is illustrated 
in Figure 5. The analytical solutions confirm the results of the qualitative analy- 
sis: the LB diffraction enlargement is significantly smaller than the enlargement 
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Figure 3: View of the spot, x, y plane, of the same long Gaussian pulse as in 
Fig. 2. The transverse size grows approximately twice over the same normalized 
distance z = Air/ 15 as in the case of an optical beam. 



of long optical pulses and optical beams. Let us remark once again that this 
result is only correct if the number of harmonics under the pulse, multiplied by 
27r, (dimensionless parameter a) is large. 

7.3 Dynamics of of light disks in linear regime. Difrac- 
tionless pulses 

As we mentioned in the beginning, producing optical pulses with small longitu- 
dinal and large transverse size, while at the same time keeping a large number 
of harmonics under the pulse, no longer presents significant experimental dif- 
ficulties. This can easily be realized in the optical region for pulses with time 
duration from 200-300 fs up to 50-60 fs. We will again consider the propagation 
of LD in the framework of the solutions l|71|) of the VLAE (|47|1 in a Galilean 
coordinates under initial conditions in the form: 



A = A x x; a = 40; S 2 = = 81, 

A x (x, y, z',t' = Q) = exp (- + £ + - ) . (82) 

The results of the calculations of the solutions (|71|l using FFT and inverse FFT 
are presented in Figure 6. Again, as one should expect bearing in mind the 
analysis in the previous section, the diffraction widening is practically absent. 
The evolution of the optical disk presented in Figure 6. demonstrates that its 
shape is preserved over the normalized distance and for time t = 375tt. This 
is larger by four orders of magnitude, than in the cases of diffraction of light 
beams and long optical pulses. This is why we refer to optical disks as being 
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Figure 4: Evolution of a Gaussian LB under initial condition A x (x, y, z,t = 0) = 
cxp (- ar+ 2 2+z ~ ) , a = 40; 5 = 1 governed by the VLAE JUJ). We present the 
plane of the pulse. A significantly smaller diffraction enlargement is observed, 
as compared with the cases of LF and light beams. 

diffractionless pulses. The numerical experiments demonstrate that such pulses 
can propagate in media with week dispersion, in dispersionless media and in 
vacuum while conserving their shape over a distance of more than hundred km. 

8 Conclusion 

The method applied of amplitude envelopes give us the possibility to investigate 
and compare the propagation of optical pulses in media with strong dispersion 
with this in media with week dispersion, in dispersionless media and in vac- 
uum. In the case of media with dispersion, we obtained an integro - differential 
nonlinear equation, describing propagation of optical pulses whose time dura- 
tion is of the order of optical period of the carrier frequency, and also of pulses 
with many harmonics under the pulse. In the case of slowly varying ampli- 
tudes (many harmonics under the pulse) we reduced this amplitude integro - 
differential equation to amplitude vector nonlinear differential equations and 
obtained different orders of dispersion of the linear and nonlinear susceptibil- 
ity. In the second case, propagation of optical pulses in dispersionless media 
and vacuum, we obtained an amplitude equation witch is valid in both cases, 
namely, pulses with many harmonics and pulses with only one-two harmonics 
under the envelope. We normalized these amplitude equations and obtained five 
dimensionless parameters determining different linear and nonlinear regimes of 
propagation of the optical localized waves. For optical pulses with small trans- 
verse and large longitudinal size (optical filaments) we obtained the well-known 
paraxial approximation, while for the case of optical pulses with a relatively 
equal transverse and longitudinal size (the so-called light bullets) and for the 
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case of large transverse and small longitudinal size (light disks), we obtained 
new non-paraxial nonlinear amplitude equations. When the optical field is with 
low intensity, we reduced the nonlinear amplitude vector equation governing 
the evolution of LB and LD to linear amplitude equations. Surprisingly, in lin- 
ear regime the normalized amplitude equations for media with dispersion and 
the amplitude equations in dispersionless media and vacuum are identical with 
accuracy of material constants. The linear equations were solved using the of 
Fourier transform technique. One unexpected new result is the relative stability 
of LB and LD and the significant reduction of the LB and LD diffraction en- 
largement with respect to the case of long pulses in linear regime of propagation. 
It is important to emphasize particularly the case of LD which are practically 
diffractionlcss over long distances exceeding hundred kilometers. 
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Figure 5: Evolution of the same Gaussian LB as in Fig. 4 , governed now by 
the VLAE in a Galilean coordinate frame (|47[) . The xy plane of the pulse 
is presented. The LB doubles its transverse (and longitudinal) size over the 
normalized distance-time t' = z' = 15tt. This is approximately 10 2 larger than 
the diffraction length of LF and light beams. 
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Figure 6: Evolution of an optical disk governed now by the VLAE in a Galilean 
coordinate frame H47JI . The xy plane of the pulse presented. The LD preserves 
its shape over a distance-time 375-7T, which is larger by 10 4 than the standard 
diffraction length Zdiff — k$r\ of light beams and long optical pulses. We call 
this type of optical pulses diffractionless. 
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